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* 20.  (Continued) 


> band  shape  of  the  measured  band  with  theoretical  spectra  syn- 
thesized at  different,  but  specific,  rotational  temperatures. 

Three  methods  of  comparison  were  used  and  it  was  concluded  that 
a sum  of  s()uares  weighting  approach  is  near  optimal  for  estima- 
ting rotational  temperatures . The  effect  of  noise  was  tested  by 
adding  real  typical  noise  to  the  synthetic  model.  The  same  com- 
parison technique  was  then  applied  to  assess  the  benefits  of 
noise.  The  uncertainty  of  the  measurement  is  estimated  to  be 
+ 5'’K.  A 'wavelike*  fluctuation  of  the  apparent  rotational  temp- 
eratufc,  with  a period  on  the  order  of  1/2  to  1 hour,  was 
observed.  The  fluctuation  on  this  occasion  appeared  to  covary 
with  the  intensity  of  the  OH  airglow  emission.  | 


LIST  OF  COnTRlBUTORS 


Doran  J.  Baker  - Principal  Investigator 

U.  Ralph  F.nibry 
Will iain  R.  Pendleton 
Allan  J.  Steed 
Gene  A.  Ware 


ACKNOWLEDGMLNTS 

The  success  of  this  experimental  program  was  made  possible  through 
the  combined  efforts  of  numerous  membf'rs  of  the  staffs  of  the  Defense 
Nuclear  Agency  (program  sponsors),  the  Air  Force  Geophysics  Laboratory, 
the  University  of  Alaska,  Utah  State  University,  and  many  others 
part  icipatc'd  in  the  ICFCAP  '7b  measurf'ments  program  at  Poker  Flat.  Al.isk,). 
liu’se  contributions  are  gratefully  acknowledged.  The  efforts  of  Dr.  A.l. 
Stair,  Jr.,  arul  Mr.  J.G.  illwick  (both  of  AIGI  ),  Dr.  k.D.  Baker,  Mr.  A.J. 
Steed.  Mr.  G.A.  Ware.  Mr.  II. R.  Lmbry  and  Dr.  W.R.  Pendleton  (all  of  USD) 
on  behalf  of  the  pro.ject  are  specifically  acknowledged. 


I ) I 


TABLE  OF  CONTENTS 


Page 

LIST  OF  CONTRIBUTORS  AND  ACKNOWLEDGMENTS  iii 

RELATED  CONTRACTS  iii 

LIST  OF  TABLES vii 

LIST  OF  FIGURES ix 

INTRODUCTION  1 

Background  1 

Emission  Mechanisms  6 

OBJECTIVES 11 

SYNTHETIC  SPECTRUM  13 

Line  Spectrum 13 

Synthetic  Spectrum  14 

INSTRUMENT  DESCRIPTION  21 

SPECTRAL  COMPARISON  METHOD  25 

Comparison  Technique  25 

Recommendations 30 

CONCLUSIONS  AND  RECOMMENDATIONS  39 

Conclusions 39 

Recommendations  39 

REFERENCES 41 

APPENDIX  A - Notation  of  Energy  State  of  Molecule  A-1 

APPENDIX  B - Relative  Abundance  of  0 (A'^E"  v=l), 
with  Respect  to  Temperature.  B-1 

APPENDIX  C - Theoretical  Review  C-1 

Vibrating  Rotator  Model  C-1 

V 


KHECKDIi'G  PAjK  tLANK-UOT  /IL-I-JD 


TABLE  OF  CONTENTS  (cont.) 


Hund's  Coupling  

Transitions  of  Oalfc'r^  ->■  

APPENDIX  0 - List  of  Computer  Programs  

APPENDIX  E - Values  of  Instrumental  Response  

APPENDIX  F - Measured  Zenith  Emission  Rates  of 

Atmospheric  (0,1)  Band  

APPENDIX  G - Measured  Spectrum  of  02(b*i:i)  Atmospheric 
(0,1)  Band  at  Poker  Flat  on  March  1,  ly75  

APPENDIX  H - Distribution  List 


C-1 

C-5 

D-1 

E-1 

F-1 

G-1 

H-1 


LIST  OF  TABLES 

Table  Page 

1.  Daytime  atmospheric  airglow  species  2 

2.  Band  systems  and  wavelengths  of  transitions 

of  O2 3 

3.  The  band  systems  of  0;  in  the  airglow  spectrum 4 

4.  O2  atmospheric  band  radiances  under 

auroral  conditions  ^ 

5.  Rotational  temperatures  of  02{b'-Zg)  (0,1) 

band  calculated  using  three  approaches  32 

6.  Ratational  temperatures  computed  from  02ib^Z^) 

atmospheric  (0,1)  band  at  8645  A . . . 36 

C.l.  Ratational  constants  of  the  X^Zz  and  b^Zg 

states  of  O2 

C.2.  Wavenumbers  (cm"*)  of  lines  in  (0-1)  band 

of  02(^Z  - *Z) C-9 


1 

vi  i 


1 


LIST  OF  FIGURES 


Figure 

1.  Simplified  energy  levels  and  transitions 

of  0? 

2.  Typical  average  altitude  profiles  for  [Oj], 

[0i]  and  temperature  for  winter  conditions 
at  60'’  latitude 

3.  atmospheric  (0,1)  band  synthetic 

line  spectrum  at  T = ZIO^K  

4.  The  synthetic  spectrum  of  Ozib^T.*)  atmospheric 

(0,1)  band  with  3.0  cnr'  resolution  

5.  The  synthetic  spectrum  of  02(/;‘S^)  atmospheric 

(0,1)  band  with  4.34  cm"'  resolution  

6.  Theoretical  spectra  of  02(Tj*F+)  atmospheric 

(0,1)  band  at  different  tetfiperatures 

7.  A simplified  schematic  diagram  of  the  EDL 

interferometer-spectrometer  system  

8.  The  night  airglow  spectrum  of  Ozib'):*]  atmospheric 

(0,1)  band  measured  at  Poker  Flat,  Alaska  at 
8:45  hrs.  UT(x=121°)  on  March  1,  1975  

9.  The  relative  intensity  of  sunlight  ana  the 

instrument  response  with  respect  to 

sunlight  

10.  The  normalized  total  response  with  respect  to 

wavenumber  over  the  range  of  02(i<‘2^)  at- 
mospheric (0,1 ) band 

11.  Comparison  of  theoretical  models  of  different 

temperatures  and  measured  spectrum  

12.  Apparent  rotational  temperature  of  0?(7)‘z^) 

atmospheric  (0,1)  band  of  each  scan  .'..... 

13.  Theoretical  volume  emission  rates  for  the 

02(/‘,.+  ) atmospheric  (0,1)  band  night 
ainjlow 

14.  Day  ainilow  volume  emission  rates  profile  of  the 

0 (.  '>’*;)  atmospheric  (0,1)  bana  from  rocket 
borne  photometer  


Page 

5 

8 

15 

17 

18 
19 
22 


26 

27 


29 

31 

33 

35 

35 


1 X 


rTicICiiHI.O  I’A";-.  -LANK-IiOT  e’lLI-illD 


LIST  01-  FIGURLS  (cont.  ) 


1 ujure  Page 

15.  Apparent  correlation  between  the  rotational 
temperature  of  0;> (; 'i;!)  atmospheric  (0,1) 
band  and  the  zenith  radiance  of  the  0H(5,3) 
band  observed  above  Poker  Flat,  Alaska  on 

March  1 , 1975  37 

C.l.  The  vector  diagram  of  Hund's  case  (;  ) C-4 

C.2.  Transition  between  and  ,v'F" C-8 

F.l.  The  observation  of  zenith  emission  of 

02(.''‘):f)  atmospheric  (0,1)  band  from  rocket F-1 

F. 2.  Zenith  emission  rates  observed  by  rocket 

borne  photometer  and  spectrometer  F-1 

G. l.  Coadded  ('.')  spectrum  0.- (;  '?:t)  atmospheric 

(0,1)  band  obtained  from'zenith  observations 
at  Poker  Flat,  Alaska  at  0831-0834  hrs.  UT, 

on  March  1 , 1975  G-1 


16 


1 


INiRODUCTION 


Background 

Many  atomic  and  molecular  species  contribute  to  the  radiation  from 
the  upper  atmosphere.  A partial  list  of  the  observed  day  airglow 
radiations,  made  at  high  altitude,  is  shown  in  Table  1.  Most  of  these 
cannot  be  observed  at  the  surface  of  the  earth,  because  the  intensity 
of  scattered  sunlight  is  much  greater. 

During  twilight  and  nighttime,  surface  observation  of  the  emission 
from  several  upper-atmospheric  airglow  species  is  possible.  Of  par- 
ticular interest  here  is  the  radiation  from  excited  oxygen  molecules. 

The  atmospheric  band  system  of  0?,  was  first  discovered  by 

:-‘eincl  [1950,1'JSl]  in  the  airglow  and  later  in  aurora.  (For  molecular 
state  symbolism  see  Appendix  A). 

The  spectra  of  atmospheric  and  infrared  atmospheric  bands  of  Oj  are 
now  well  established  as  features  of  the  airglow  spectrum.  A deta 
listing  of  the  band  systems  for  known  transitions  of  O2  is  given  in 
Table  2.  Table  3 gives  the  observed  radiances  for  the  Herzberg  bands 
and  two  bands  of  the  atmospheric  and  infrared  atmospheric  band  systems, 
all  taken  at  high  altitude  using  bal loon-borne  eguipment. 

Figure  1 shows  a simplified  energy  level  diagram  for  the  oxygen 
molecule,  with  two  vibrational  levels  corresponding  to  v=0,  and  v=l . 

The  (0,0)  band  at  7600  A and  the  (0,1)  band  at  8645  A correspond  to  the 
transitions  from  the  state  (.'*;■*■  , v=0)  to  the  states  (v’il',  v=0)  and 
(A'^'I*’,,  v = 1),  respectively.  According  to  the  airglow  radiances  of  Table 
3,  and  the  auroral  radiances  given  in  Table  4,  the  radiance  of  the  (0,0) 
band  is  typically  about  20  times  that  of  the  (0,1)  band.  In  ground- 
based  observations , the  (0,0)  band  is  completely  reabsorbed  by  the 
great  amount  of  0.  in  the  ground  state  (a’2',  v=0)  in  the  lower  atmos- 
phere. The  low  concentration  of  O.  molecules  in  the  vibrational ly- 
excited  state  v=l)  causes  little  absorption  of  the  (0,1)  band. 


..vTOBB.'vrr’v  . 


17 


Table  1.  Daytime  atnwspheric  airglow  species  \_Whittcn  and  Poipofj',  19 


Emission 

Wavelength 
(A  or  uni) 

Zenith  Radiance 
(kR) 

Lyman  o 

1216 

5-12 

01  ( 'P- 

1304 

2-6 

01 

1355 

0.4 

N0(v) 

2000-3000 

1 

N (2PG) 

3000-4000 

0.4 

Nt(lN) 

3914 

2-7 

Ni(\;-p-) 

5200 

0.1 

0I(‘.  -'  ■) 

5577 

0.4-3 

Na(.  ) 

5893 

2-40 

0I(^^-‘  ) 

6300 

3-60 

o,(/  ’ t-A- 

7600 

'300 

8640 

~ 20 

1.27 

'30,000 

OH 

2. 8-4.0 

'5,000 

I 


able  2.  Band  systems  and  wavelengths  of  transitions  of  0,  and  1962]. 


' Table  3.  The  band  systems  of  O2  in  the  airglow  spectrum  \_Evans  and 

Llcwp.lJyn,  1970]. 

I 


Day  Airglow 

Night  Airglow 

i 

I Band  System 

Transition 

Wavelength 

Radiance 

Radiance 

(A  or  pm) 

(kR) 

(kR) 

1 Herzberg  bands 

- X^l- 

3100  - 5000 

No  observation 

0.5 

u g 

. 

i Atmospheric 

- X^Y.- 

7619  (0,0) 

300 

8 

1 System 

'J  9 

8645  (0,1) 

15 

0.4 

Infrared 

a'A  - A'^Z- 

1.27  (0,0) 

( 

{ 

\ 

Atmospheric 

System 

9 9 

1.58  (0,1) 

25000 

90 

Table  4.  O2  atmospheric  band  radiances  (kR)  under  auroral  conditions, 
1 (5577)  = 100  kR  [Vallanec  -forif-s  and  Jatt'ngrv,  1974], 


y 0 1 2 3 4 5 


0 

7619 

8645 

(1200)+ 

(58)* 

1 

6882 

7708 

8742 

(13)* 

(1.6) 

2 

6286 

6968 

7802 

8844 

(1.9)* 

(6.8)* 

(1.4) 

3 

6370 

7059 

7901 

8950 

(0.40) 

(6.2)* 

(12) 

(3.9) 

4 

645/ 

7154 

8006 

(0.30) 

(2.8)* 

(3.1) 

5 

6457 

7255 

8118 

(0.55) 

(3.5) 

(2.3) 

he  wavelength  for 

each  band  is  (jiven  in 

angstroms . 

The  value 

in  paren 

theses  are  corresponding  radiance  values. 


*L)ncerta i nty  less  than  25  . 

+Derived  from  1(0,1)  from  transition  probability. 


A 


b 


By  using  the  Boltzmann  probability  distribution  formula,  the  ratio 
of  abundance  of  0*(X^T.~,  v=l ) to  v=0)  can  be  well  approximated 
by  the  expression: 


where 

K = Planck’s  constant,  6.625  x 10"^^  erg-sec 

= Velocity  of  light  in  vacuum,  2.99793  x 10'°  cm/sec 

\>i  = Wavenumber  of  the  transition  from  the  first  vibrational 

state  (A'^2",  v=l ) to  the  ground  state  (.v^2",  v=0),  1556  ciir ' 

k = Boltzmann's  constant,  1.38044  x 10"“'  erg/deg 

T = Temperature,  in  °K 

Values  of  r(.)  for  various  temperatures  are  listed  in  a table  in  Appendix 
B.  In  the  tropospheric  condition,  r(7')  is  always  less  than  about  2 x 10" 
Although  some  measurements  have  been  taken  by  rocket  observations 
and  1968],  there  is  considerable  advantage  in  being 

able  to  economically  measure  the  emissions  of  0.  in  the  mesosphere  from 
ground-ba  ed  observatories  over  a long  period  of  observation. 

Ihis  study  was  done  on  the  (0,1)  band  of  the  atmospheric  system  of 
0.  . The  data  analyzed  here  were  taken  jy  USU  personnel  in  March,  1975 
using  the  USU  Arju.-  Mobile  Aeronomy  Observatory  which  at  the  time  was 
located  at  Poker  Flat,  Alaska. 

Enns^i 0 n Mecha nisms 

In  their  paper  dealing  with  rocket  observations  of  the  atmospheric 
bands  in  the  day  airqlow,  Wa!!a,<,  and  lIunt’K  [1968]  suggest  three  sources 
of  0..  (.•  which  would  be  expected  to  occur  only  in  a sunlit  atmosphere. 
They  are: 

1.  Energy  transfer  from  0('  ),  dominant  above  100  km 
1300-1  750  A)  + 0.  • 0(  ’P)  + 0( ) 

0('; ) 1 0,  (,v^:  -,)  • o(  ) + 0. 


(2) 

(3) 


2.  Resonance  scattering,  dominant  in  the  region  from  65  km  to 
100  km 


; -(7600  A)  ♦ 0 • o,.(;  ‘^:+) 


3.  Photoly. is  of  ozone,  dominant  in  the  region  fror,  35  km  to 
65  km  and  probahly  imiiortant  u[-  to  100 


. -(  - A 2300-3000  A)  • 0 . • 0 ( 1 


0(‘  ) ^ 0, (.  -) 


) * o.(  • V 


Later  work  has  shown  that  : he  r-u'-tu-i'  den-  t-  of  excited  oxygen 
molecules  is  closely  related  to  tiu-  rumiber  den-  yy  uf  ozon«.  Reaction 
(4)  can  be  repl.iced  by  [ • ' and  • , l'-*72i 

• ('■  2000-26O0  a;  * U • 0 (-'■*)’  0i  (5b; 

2660-3300  A)  f 0.  • 0.  ( • , ' 0(  \ ) (5.  1 


which  seems  to  dominate  produiition  in  daylight  from  40  to  100  km 
altitude.  Figure  3 stiows  .iltitude  profiles  for  0,,  0.  and  the 
temperature  for  arctic  winter  conditions. 

Botii  morning  and  evening  conditions  can  be  considered  as  a trans- 
ition period  between  the  day  airglow,  in  which  resonance  radiations  are 
dominant,  and  the  night  airglow.  in  which  chemiluminescence  becomes 
important.  The  morning  twilight  enhancement  iireatly  exceeds  that  m 
the  eveninc)  in  the  80  to  100-km  region  [ '.  . - ■ . 1')7‘jJ;  (he  reason  for 
this  is  the  (iroduction  of  0.  (.  M l>y  energy  (r.insf'er  from  0(  • ) 
produced  in  ozone  photolysis  at  the  '.ame  altitude.  At  sunset,  the 
production  of  0(  ’ ) in  0^  ptiofoly.is  occurs  only  above  lt>0  km  [ '.  ■ 

and  , 1(721,  and  tti'-  ri’sultim]  yield  of  0 (-''^)  is  negligible. 

In  ♦ he  night  airglow,  the  absence  of  sunlight  renders  chimiical 
and  ion-electron  reeomb  i na  t i on  r'eactions  ttu-  don-inant  souri  es  oi 
0 (.-  ‘ ^).  The  resonance  scattei’ing  that  predominate-  during  daytime 
and  tw'i  light  conditions  is  relegated  to  ,i  minor  role.  the  excited 
states  ot  moleculiit-  oxygen  i an  be  forrneil  by  ttuee-boily  recomtu  nat  ion ; 


CONCENTRATION  (cm' 


f 


i 

\ 

1 

i 

I 


ro 


2 

UJ 

CO 

c 

o 

1- 

o 

<TJ 

o 

<V 

lO 

C 

fV 

o 

o 

0) 

> 

S- 

a; 

(UJ>()  1H9I3H 


9 


o(  '/■)  + o(  '/■)  ^ M > o..(/  ■;;+)  + M (6) 

Under  auroral  conditions  the  iiietastable  states  of  molecular 
oxygen,  0.  (-  in  particular,  can  be  produced  by  secondary  auroral 

electrons,  that  is; 

. - + 0. ► . ' + QiU-'y*)  (7) 

'J 

The  0r(.  ‘:  '*■)  can  also  be  formed  under  auroral  conditions  by  thermal 
col  1 isions  wi th  metastable  oxygen  atoms: 

o('; ) + 0? • 0(  '/>)  + 02(/  'i:^^)  (3) 
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OBJECTIVES 

The  objective  of  this  study  is  to  estinidte  the  rotational 
temperatures  of  0?(/ ► A'’’):")  atmospheric  (0,1)  band  at  8645  A 
in  the  airglow  by  the  following  steps: 

1.  Obtain  the  observed  spectra  of  the  night  airglow. 

2.  Compute  the  synthetic  spectra  at  different  temperatures 
for  an  appropriate  theoretical  model. 

3.  Optimize  a technique  for  comparing  observed  spectra  with 
synthetic  spectra  of  different  rotational  temperatures. 

4.  By  the  comparison  of  the  observed  spectra  with  synthetic 
spectra  ascertain  the  apparent  rotational  temperatures 
of  0p(''C^)  radiation  of  the  night  airglow. 


tLAlIK-ncT  c-lL-siD 
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SYNTHETIC  SPECTRA 

U ne  Spectrum 

The  radiance  of  an  emission  line  in  the  infrared  rotational 
spectrum  or  rotation-vibrational  spectrum  not  only  depends  on  the 
Boltzmann  distribution  function  but  also  depends  on  the  quantum 
number  of  the  upper  states.  The  emission  radiance  for  a given  line 
can  be  expressed  as:  \_Herzbei?g , 1950] 


2.'  r + 

r - I'-W  „ U' 

I r X 6 . , X exp 

em  Q ’!  ^ kT 

i‘  L 


where 


= The  radiance  of  an  emission  line,  in  photon  sec'‘cm'^sr"‘ 
= The  emission  constant 

■j  = The  wavenumber  of  each  transition  which  has  been 
calculated  and  listed  in  Table  C.2 

■;  = The  partition  function 


' (2-+1)  X exp  [-“'(•')  X Cc'/P.’’] 


■ ' = The  total  quantum  number  of  the  upper  state 
.S'  , = The  line  strength  of  each  transition. 

ej 


The  ratio  of  splittino  constant  A to  the  molecular  rotational 
constant  nameiv.  a I 4 Appendix  C).  _ , tan 
be  exiiressed  as  follows  , 1936]: 


for  branch 


for  ^0  branch 
for  branch 


for  branch 


rHKCliDIi'G  i^AGE  bLANK-MOT  fILi-1.!:D 


Tills  fiildtionship  is  based  upon  an  assumption  of  rotational  equilibrium 
in  the  upper  state.  One  should  justify  this  for  the  system  of  interest. 

In  this  paper,  only  the  relative  radiance  is  needed.  The  attempt 
is  by  comfiarison  of  the  shapes  of  theoretical  band  profiles  with  measured 
profiles  to  deduce  the  apparent  rotational  temperature.  The  emission 
constant  ' „ in  Equation  (8)  need  not  be  evaluated  in  the  calculation 
of  relative  radiance.  The  relative  radiances  of  line  spectrum  at  a 
temperature  of  T=210"K  are  shown  in  Tiqure  3.  This  band  shows  two 
winqs,  the  P wing  which  includes  branch  and  branches  and  an  R 
wiruj  which  includes  an  branch  and  a branch.  The  band  center  is 
at  11564  cm" ‘ . The  computer  programs,  named  02LMI  and  SORT,  which 
generate  the  wavenumber  and  radiance  of  each  line  and  sorted  in  the  order 
of  increasing  wavenumber  respectively,  are  listed  in  Appendix  D. 

Spectrum 

For  each  particular  wavenumber  point  of  the  band,  the  radia  ce 
can  be  computed  by  the  formula  as: 

sire  . . 


where 

Ay..  = •e.-'j.|,  in  cm"  ‘ 4xRES 

^ Ihe  wavenumber  of  any  point  ' in  the  hand  in  cr  * 

. = The  wavenumber  of  one  line  in  line  spectrum  in  cm"' 

fhe  resolution,  which  means  the  wavenumber  interval 

from  the  central  maximum  to  first  minimum  of 

s i n ■ .. 

function 

- Thi>  rel.ifive  radiance  of  one  lino  in  the  line  spectrum 

X . 

, b 

Total  number  of  emission  lines,  namely,  77  for  . 0 to 
RLS  - Resol ut ion 


i-r-ifiCiiDI.G  I’Ajt  ELA!lK-noT  i-'ILi-i.iD 


synthetic  line  spectrum  at 
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The  apodizinq  function  = 1 


2 


Taking  into  account  the  resolution  function  and  the  apodizinq  effect 
of  the  filter,  it  is  assumed  that  the  value  of  the  spectral  radiance  can 
be  neglected  beyond  the  fourth  minimum  on  each  side  of  the  function 

contributed  by  the  radiance  of  a single  emission  line.  The  spectral 
radiance  at  each  wavenumber  point  of  the  band  region  is  the  sum  of  all  the 
contributions  at  this  point.  The  relative  radiance  of  the  whole  band  has 
been  calculated  with  a 0.5  cm'*  interval  over  the  band  region  for  com- 
parison with  observed  spectra. 

The  computer  program,  named  CNVLS,  is  listed  in  Appendix  D.  The 
spectral  radiance  profiles  for  different  temperatures  and  resolutions 
are  shown  in  Figures  4,  5,  and  6.  In  Figures  4 and  5,  the  envelope  of 
the  line  spectrum  becomes  smoother  as  the  resolution  interval  increases. 
Two  spectra  for  T=150°K  and  T=250"K,  both  at  3.0  cnri  resolution,  are 
shown  in  Figure  6.  The  peak  of  R wing  and  the  peak  of  P wing  each  shift 
away  from  the  band  center  when  the  temperature  increases. 


Figure  4.  The  s/nthetic  spectrum  of  0,  ( ')  atmosphr-'' ic  (0  J ) band  with  3.0  cnr  • resolution. 


WAVENUMBER  (cm* 
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INSTRUMENT  DESCRIPTION 


The  measuring  s>"tem  for  the  airglow  emission  is  a large-aperture 
field-widened  interferometer-spectrometer  system  which  has  been  designed 
and  developed  by  the  Electro-Dynamics  Laboratories  of  the  Utah  State 
University  \_Deouain  el  al. , 1970].  A simplified  schematic  diagram  of 
this  interferometer-spectroihoter  system  is  shown  in  Figure  7.  Two 
wedge  end  mirrors  are  used  instead  of  the  simple  mirrors  of  a conven- 
tional Michelson  interferometer.  One  of  the  wedge  end  mirrors  is  fixed 
and  the  other  is  driven  by  a mirror  carrier  system  which  includes  a 
laser  beam  interferometer  to  record  the  position  of  the  movable  wedge 
end  mirror.  The  salient  features  of  the  system  are:  (1)  field  of 

view  widened  to  5°  full  angle;  (2)  a 5-cm  clear  aperture;  and  (3)  a 
5-cm  drive  capability. 

Suppose  that  a beam  of  light  is  divided  by  the  beam  splitter  into 
two  mutually  coherent  beams,  as  in  the  Michelson  interferometer,  and 
these  beams  are  reunited  after  traversing  different  optical  paths.  The 
reunited  beam  goes  through  a lens  and  is  detected  by  a photodetector. 

If  the  light  is  not  monochromatic  but  has  a spectral  composition  given 
by  the  power  spectrum  <;(v),  where  u is  the  wavenumber,  then  the  in- 
tensity of  the  reunited  beam  varies  in  a manner  that  depends  on  the 
particular  optical  spectrum.  By  recording  the  intensity  as  a function 
of  the  path  difference,  or  the  position  of  the  movable  wedge  end  mirror, 
the  power  spectrum  ;(u)  can  be  deduced.  We  can  write  the  intensity  of 
the  reunited  beam  ;(x)  as: 


or 


(11) 


.V(.r-)  ^ 2l{x) 


(12) 


wLAIJK^IJOT 
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where  /(O)  is  the  intensity  for  zero  path  difference  and  x is  the  path 
difference.  Therefore,  iv'(x)  and  j{j)  constitute  a Fourier  transform 
pair.  Accordingly,  we  can  write; 


;(v)  = y v(x)  dx-  (13) 

-oo 

From  the  output  of  the  detector,  a(.v)  is  the  output  voltage  signal 
which  depends  on  the  position  of  the  movable  wedge  end  mirror  recorded 
by  the  laser  as  a reference  signal,  [hen,  the  Equation  (13)  can  be 
solved  by  the  Fourier  transform  computer  program. 

The  resolution  capability  of  a Michelson  interferometer  is  gener- 
ally taken  to  be; 


Av  = ^ 


where  Av  = vp  - vj,  is  the  difference  in  wavenumber  between  two  re- 
solvable adjacent  monochromatic  lines  and  6 is  the  maximum  optical  patti 
difference.  The  relationship  between  the  optical  path  difference  and 
distance  of  motion  of  the  wedge  prism  along  its  apparent  mirror  image 
plane  in  the  field-widened  system  can  be  expressed  as  1 976]: 


= 2x  }i  sin  a + tan  y 


>1  ^ s i n " g - 1 
cos  b 


cos  y 


where  .r  is  the  drive  length,  n is  the  index  of  refraction,  u is  the 
prism  angle,  p = sin"‘(H  sin  u) , and  y = tan~'  (tan^a/tan  p).  For  the 
values  u = 8°  and  = 1.5; 

cS  = 0.2b-  (16) 

The  drive  length  x can  be  derived  as; 


= ,v  ^ X y - X , 
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where  .V  is  the  number  of  cycles  of  the  reference  signal  between  each 

L? 

sample,  is  the  wavelength  of  laser  and  w is  the  total  number  of 

L ^ 

sample  points  in  an  interferogram.  For  the  values  .V  _=  4,  X^=6328  A and 

-V  =16384,  the  drive  length  u;  is  equal  to  2.0735  cm  or  6=2.296  cnr*. 
s 
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SPECTRAL  COMPARISON  METHOD 

Comparison  T echnique 

While  recording  the  interferograms,  the  incoming  sky  radiation 
was  limited  on  the  long  wavelength  side  by  the  response  of  the  photo- 
detector, and  on  the  short  wavelength  side  by  a filter.  Figure  8 shows 
a typical  spectrum  obtained  by  digitizing  the  interferogram,  computing 
the  spectrum  by  using  the  fast  Fourier  transform  (FFT)  and  plotting  the 
amplitude. 

Before  using  such  a spectrum,  it  must  be  corrected  for  the  instrument 
response.  This  was  accomplished  by  obtaining  the  spectrum  of  a ra- 
diation source  of  known  characteristics,  and  ascertaining  a correction 
factor  as  a function  of  wavelength. 

An  interferogram  was  recorded  with  the  instrument  pointing  at 
the  full  moon,  and  the  corresponding  spectrum  was  computed  and  plotted. 

The  source  spectrum  was  treated  as  a blackbody  at  a temperature  of 
6000°K,  which  the  solar  spectrum  closely  approximates  in  this  wave- 
length region.  From  Planck's  blackbody  distribution  equation: 


1 


Q'n'na  1 

^ expl^-TW^  - T 


(18) 


where 


= Energy  density  per  unit  volume  per  unit  wavelength,  in  erg/cm'*. 


The  relative  intensities  for  different  wavelengths  of  solar  radiation 
can  be  derived  by: 


1 ) 


A , 


X 

X 


exp  {lu-/k\-y:)  - 1 
exp  (l;c*/A-A ) '/■)  - 1 


(19) 


figure  9 shows  a smooth  plot  of  a portion  (11457  to  11614  cm"') 
the  spectrum  obtained  from  observing  the  full  moon  with  the  same 


spectrum  of  atmospheric  (0,1)  band  measured  at  Poker 


Figure  '5.  Fhp  rplativp  intonsit.  nt  s.jntigh^  and  instruncnt  responsp  with  respect 
to  sjnl ighr  . 
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spectrometer.  Also  shows  a plot  of  the  relative  spectral  radiance  of  sun- 
light modeled  as  a 6000°K  blackbody  source.  The  instrument  correction  for 
any  wavenumber  in  this  region  is  the  ratio  of  the  relative  spectral  radiance 
of  the  source  to  the  relative  intensity  of  the  observed  spectrum,  which  is 
shown  in  Figure  10.  The  values  of  instrument  response  corresponding  to 
each  sample  point  over  the  band  are  listed  in  Appendix  E. 

The  spectral  radiance  at  each  sample  point  of  the  spectrum 

is  instrument  corrected  by  dividing  the  instrument  response  (Fig.  10)  and 
then  is  normalized  by  dividing  by  the  total  radiance  of  the  band.  For  the 
spectrum  given  in  Figure  8 for  8:45  hrs.  LIT  (solar  zenith  angle  x = 211  °) 
taken  at  Poker  Flat,  Alaska,  on  March  1,  1975,  the  apparent  zenith  radiance 
of  the  02(7)^':''')  (0,1)  airglow  band  was  500  rayleighs.  This  value  is 
uncorrected  for  atmospheric  extinction.  Then  each  corrected  sample  value 
is  treated  as  radiance  of  one  line  in  the  theoretical  line  spectrum,  and 
interpolated  by  — function  with  a resolution  equal  to  3.0  cm"'. 

A family  of  theoretical  spectra  has  been  generated  over  the  range 
from  140  to  300°K  in  10°K  steps.  These  were  convolved  to  a resolution  of 
3.0  cm"'.  Values  of  the  spectral  radiance  of  both  observed  and  theoretical 
spectra  were  calculated  at  sample  points  spaced  each  0.5  cm* ‘ over  the 
range  from  11400  cm"'  to  11650  cm"'. 

The  range  of  wavenumbers  for  optimal  band-shape  comparison  were  limited 
to  11535.9  to  11616.3  cm"'  because  of  overlap  with  other  spectral  features 
of  the  airglow  and  aurora  [WalLat^e  and  Chtvv;  rrLiLn,  1 959].  These  include 
multiplet  of  NI(11431 .2-11517.3  cm"'),  the  (2,1)  band  of  the  First  Positive 
system  of  N?  in  the  region  11461.8-11530.6,  and  the  Pi(7)  and  Pp(7)  lines 
of  the  (6-2)  band  of  0H(11628.5  and  11640.7  cm-'). 

The  best-fit  band  shape  comparison  for  rotational  temperature  esti- 
mating proceeds  as  follows.  The  total  band  radiances  of  both  the  synthetic 
and  the  experimental  spectra  within  the  selected  comparison  range  are 
forced  to  equality  by  change  of  scale  factor.  The  difference  between  the 
radiances  of  the  observed  and  the  theoretical  spectra  were  calculated  at 
each  0.5  cnr ' in  the  range.  The  sum  of  the  absolute  values  of  the 
differences,  the  sum  of  squares  of  the  values  of  the  differences, 

and  the  sum  of  relative  absolute  values  of  the  difference, 
were  all  calculated  for  the  specified  rotational  terii()era  tures . These  vlaues 


atmospheric  (0,1)  band. 
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change  from  one  rotational  temperature  to  another,  and  the  minimum  of 
difference  was  found  by  a curve  fitting  technique  as  shown  in  Figure  11. 

The  observed  spectrum  processini]  program,  DTPRS,  the  comparison 
program  and  the  numerical  results  associated  with  Figure  11  are  listed  in 
Appendix  D.  As  expected,  the  variance  of  temperature  given  by  the  different 
comparison  techniques  was  dependent  upon  the  signal  to  noise  ratio.  The 
sum  of  the  squares  technique  was  adopted  as  optimal  because  it  empirically 
gives  the  best  weighting  of  the  various  portions  of  the  spectrum  according 
to  strength  of  signal  relative  to  noise. 

Pi scuss  ion^  o^f  Resul  ts 

All  of  the  0.(.-'..^)  spectra  considered  here  wore  taken  at  Poker  Flat. 
Alaska,  on  March  1,  1975.  Rotational  temperature  for  each  spectral  scan 
calculated  by  the  curve-fitting  technique  are  listed  in  Table  5 using  thre(> 
different  approaches  as  discussed  in  the  [)revious  section.  The  (..■)' 
approach  emphasizes  the  portions  of  the  band  which  have  higher  value  of 
spectral  radiance  in  comparison  with  the  lower  ones  which  are  more 
affected  by  the  noise,  figure  12  shows  the  rotational  temperatures  coin- 
puted  using  both  the  and  approaches.  The  results  are  generally 

in  good  agreement  with  each  other  and  also  in  comparison  with  the  result 
from  using  the  weightimj  function  Y'.‘‘ ■ where  ' is  the  radiance  of  ob- 
served sfiectrum  at  the  comparison  point.  The  agreement  is  best  when  the 
0.  {.'  is  brightest,  as  is  to  be  expected. 

The  comparison  approach  was  also  tried  on  two  coadded  (signal  averaged) 
spectra.  The  results  of  (/')'  again  show  a much  better  agreement  between 
the  apparent  rotational  temperature  of  the  coadded  spectrum  and  the 
avera<)('  rotational  tem[)erature  over  th(>  ‘.am(>  time  interval.  In  f igui'e  12, 
the  result  calculated  using  curve  shows  more  temporal  fluctuation 

than  result  from  the  other  summing  mt'tfiods.  Ifio  average  a(iparent  rotational 
temperature  over  the  time  interval  usirni  (..■■■)■  is  about  21  i K*;  this  is 
close  to  the  expected  ambient  temperature  at  about  'Ki  km  as  shown  in  ligure 

‘Points  were  iijnored  it  the  difference  between  computed  temperatures 
using  and  (.')■  was  10  K. 


Table  5.  Rotational  temperatures  ("K)  of  0.  ( ' ) (0,1)  band  calculated 

using  three  approaches. 


Time 

(hrs.  UT) 

Sum  of  Relative 
Absolute  Differences 

Sum  of  Absolute 
Di f ferences 
L|a'i  1 

Sum  Differences 
of  Square 
T.{d:.y 

8:30-8:58* 

220 

203 

221 

8:58-9:59* 

224 

201 

216 

8:31-8:34 

290 

260 

244 

8:34-8:38 

194 

193 

190 

8:46-8:50 

264 

230 

240 

8:50-8:54 

200 

192 

188 

8:54-8:58 

260 

241 

230 

9:02-9:06 

222 

230 

240 

9:14-9:18 

213 

252 

267 

9:18-9:22 

206 

199 

215 

9:27-9:31 

180 

182 

202 

9:31-9:35 

182 

180 

189 

9:35-9:39 

253 

25i 

241 

9:39-9:43 

220 

21. 

213 

9:43-9:47 

240 

235 

248 

9:47-9:51 

182 

18i 

171 

9:51-9:55 

235 

i-  J 

245 

9:55-9:59 

189 

1 

195 

8:30-8:58** 

24? 

22.) 

218 

9:02-9:59 

21 1 

21  1 

Op  j 

Over  a 1 1 Avi 

226 

21' 

*^20 

I 

1 


*Coadd  spectrum. 

**Average  temperature*  over  the  time  period. 
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Of  concern  is  whether  0?(.'  ‘),+  ) state  rotational  temperature  equilibrium 
has  been  achieved  prior  to  radiative  relaxation  by  sufficient  collisions 
with  the  ambient  atmospheric  molecules  (primarily  N2).  The  collision 
frequency  per  molecule  can  be  writtei  as; 

f = nco  (20) 


where 

n = the  density  of  air,  cm" '' 

3 'yj <5veracje  speed,  r is  Boltzmann's  constant,  is 

the  absolute  temperature,  ri  is  mass  of  molecule 
o = , the  projected  a^ea  of  a spherical  molecule; 

;■  = the  internuclear  distance,  1.904  x 10"  "cm  for  Nj, 

1.207398  X 10"  "cm  for  0?  1950] 

At  90  km,  >.  10‘‘*  cn)"\  T ~ 200'^K,  = 4.7816  x 10"‘"  g/molecule, 

and  Z‘  = 3.834  x 10"  cm/sec.  Therefore,  3.6  x 10’  /sec.  The  life-time 
of  0;>(.' '2^)  is  approx imately  to  12  sec  [.Vf.v:  1970J;  so  thermal  rotational 

j 

equilibrium  may  be  achieved. 

Theoretical  vertical  profiles  of  the  volume  emission  rate  for  the 
O'/i-' (0,1)  band  of  the  night  airglow  are  given  in  Figure  13 
1961;  .'uLi  jf' , 1966].  The  peaks  of  both  curves  occur  around 
90  km.  The  oxygen  excitation  in  the  90  km  region  is  attributed  to  the 
reaction; 


0 + 0 t M • 0,  (;  ' '.)  r M (21) 

By  comparison,  a rocketborne  measurement  of  the  daytime  airglow  volume 
emission  rate  profile  of  0,(.‘  ‘',^_)  atmospheric  (0,1)  band  as  given  by 
Wdlii-.  and  n,  [1968]  is  shown  in  Figure  14  (Appendix  F). 

Some  results  of  calculations  of  apparent  totational  temperatures  of 
the  0;(.  ';:*■)  (0,1)  atmospheric  band  by  ottu'r  workers,  as  well  as  that 
rf'ported  herein,  are  stimmari/ed  in  Tal)l(*  6. 


Og  EMISSION  RATE  (cm-=  sec"') 

Fiqurp  1.1.  Thporetical  volume  emission  rates  for 
the  02{b^Zt)  atmospheric  (0,1)  band  night  airglow. 
Curve  B is^given  by  bavth  [1961];  the  curve  labeled 
CJH  uses  the  atomic  oxygen  distribution  of 
Colegrove  et  at.  [1966J. 


VOLUMf  EMiSS.'OV  («R-  km  or  10  cm  »fc  ) 


Figure  14.  Oav  airglow  volume  emission  rates 
profile  ot  the  atmospnoric  (0,1)  band 

from  rocl'et  borne  photometer  and 
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Table  6.  Rotational  temperatures  computed  from  02(i>'^^)  atmospheric 
(0,1 ) band  at  8645°A.  ^ 


Solar 

Apparent 

Investigators 

Season  & Year  Latitude  Angle 

Condition  Temperature 

(x) 

(°K) 

Meinel  [1950] 

1949 

43° 

Quiet 

Night 

150  i 20 

Wallace  & 

Chamberlain  [1959] 

1959 

43° 

Aurora 

195  1 10 

Wallaae  & 

Chamberlain  [1959] 

1959 

43° 

Quiet 

Night 

183  * 7 

Sue  a Baker  [1976]  Spring, 

1975 

65° 

Auro  ra 

211  ± 10 

A noticeable  apparent 

correlation 

between 

the  rotational  temper- 

atures  of  band  and  the  radiance  of  the  0H(5,3)  band  was  found 

in  the  data  taken  on  March  1,  1975,  at  Poker  Flat,  Alaska.  Figure  15, 
which  shows  the  similar  pattern  of  f luctuations . These  fluctuations  might 
be  due  to  acoustical  gravity  waves  iAi’ri.'tvona,  1975],  possibly  associated 
with  the  deposition  of  energy  in  the  auroral  region.  The  change  of  inten- 
sity and  of  rotational  temperature  might  bo  due  to  the  vertical  dis- 
placement of  the  omission  layer  [h'akei-,  1976].  Further  study  of  these 
fluctuations  would  be  appropriate. 

The  strength  of  the  aurora  at  the  time  of  these  0;  (;  ’ 7''^^)  measurements 
was  ascertained  from  the  radiance  of  ttie  01(8446  A)  line'  winch  was  observed 
s imu  I t.ineous  ly . The  value  of  the  radiance  (Appendix  G)  of  this  line  was 
0.14  kR. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

Con£lu_y_ons 

1.  The  rotational  temperature  of  the  0.  atmospheric  (0,1)  band 

can  be  estimated  by  comparing  the  band  shape  of  the  observed  rota- 
tional structure  with  synthetic  spectra. 

2.  A sum  of  squares  weighting  method  to  best  fit  the  observed  spectrum 
with  the  synthetic  spectra  appears  to  be  a near  optimum  procedure. 

3.  Real  typical  noise  has  been  added  to  a synthetic  spectrum  and  used 
like  the  real  spectruni  to  analyze  the  effect  of  noise.  The  result- 
ing deviation  is  much  less  than  that  due  to  tiie  different  weiqiiting 
methods  (see  Appendix  G).  The  fluctuations  of  the  observed  tem- 
perature appears  to  be  caused  by  other  factors  in  addition  to  noise. 

4.  The  average  value  of  the  apparent  rotational  tt'  'perature  of  0 ( ' 
on  March  1,  1975  at  8:45  hrs.  UT  was  210  t.  with  the  band  radiance 
at  1/2  .J:. 

5.  A fluctuation  was  observed  which  appears  to  covary  with  the  in- 
tensity of  the  OH  airglow  emissions.  These  may  be  attributable 
to  atmospheric  wave  motions. 

Recommenda t i ons 

1.  For  more  accurate  values  of  the  spectral  line  radiances,  a curve- 
fitting technique  needs  to  be  developed  and  applied  to  the  sample 
points  near  each  line  position. 

2.  The  wavenumber  spacing  between  each  sample  point  depends  on  the 
position  where  the  wavenumber  is  located  in  the  spectrum.  An 
interpolation  function  needs  to  be  developed  which  will  improve 
the  alignment  between  the  observed  spectrum  and  the  theoretical 
spectrum. 

3.  Establish  which  part  or  parts  of  the  band  are  most  sensitive  to 
the  rotational  temperature  changes.  By  understand i ng  this,  the 
accuracy  of  the  rotational  temperature  miqht  be  able  to  tie  improved. 

4.  Take  and  analyze  data  under  both  quiet  and  known  auror.il  conditions 
and  also  at  midlatitude  to  optimize  the  model  and  minimize  effects 
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from  contaminating  species. 

5.  Make  measurements  simultaneous  with  OH  rotational  temperature  measure- 
ments. 

6.  Develop  a noise  analysis  to  better  establish  the  uncertainty  of  the 
measurement. 

7.  Establish  the  cause  of  the  apparent  "wavelike"  variations  in  the  tem- 
perature and  the  covariance  with  OH  eoiissions. 

8.  The  altitude  distribution  of  the  02('  emission  needs  to  be 
established  under  both  twilight  and  nighttime  conditions. 
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APPENDIX  A 

NOIATION  Of  ENERGY  STATE  OF  ITOLECULE 

The  notation  of  the  energy  states  of  a molecule  is  similar  to  the 
notation  used  for  the  energy  states  of  an  atom.  The  general  form  is; 


1.  A symbol  in  the  position  #1  represents  the  quantum  number  of  the 
vector  component  A along  the  internuclear  axis  of  electronic 
orbital  angular  momentum  L.  Its  magnitude  is  A(/;/20> 

A = 0,1 ,2,3, . . . ,L.  The  symbol  in  position  <>1  is; 

Y.  when  A = 0. 

II  when  A = 1 . 

A when  K - 2. 

'I'  when  A = 3. 


2.  A number  in  the  position  “2  is  equal  to  2.'  + 1 , is  the  quantum 
number  of  s|iin  S.  The  value  of  2.'  + 1 is  the  spin  multiplicity 
of  the  energy  state. 

3.  The  symbol  in  the  position  ^2  represents  the  symmetry  of  the  total 
eigenfunction  of  the  state  with  respect  to  any  plane  through  the 
internuclear  axis.  The  symmetry  is  called  "+"  if  the  eigenfunctions 
remain  unchanged  upon  coordi.iate  reflection  at  the  origin.  Other- 
wise it  is 

4.  The  symbol  in  the  position  ^4  represents  synnietry  of  the  electronic 
eigenfunctions.  It  is  it  th(>  symmetry  jiroperty  is  even. 
Otherwise  it  is  " r".  (from  Oi'rman  "gerade"  and  "un<)(n'ade" ) 

5.  Tht'  symbol  in  position  “b  indicates  the  excited  state  by 

or  The  (jround  state  is  designated  by 


p 


B-1 

APPENDIX  B 

RELATIVE  ABUNDANCE  OE  0.  ( v=l), 

WITH  RESPECT  TO  TEMPERATURE 

Relative  abundance  of  the  v=l  vibrational ly-excited  ground 
state  of  molecular  oxygen,  v=l),  as  a function  of  temperature. 


TEMPERATURE 

{''K) 


RELATIVE  ABUNDANCE 
[02]*(.t'Z-,  v=1)/[02](x'L-,  v=0) 


160 

8.382 

X 

10-" 

170 

1.909 

X 

10-'' 

180 

3.968 

X 

10-' 

190 

7.635 

X 

10-' 

200 

1.375 

X 

10-' 

210 

2.345 

X 

10-' 

220 

3.807 

X 

10-' 

230 

5.926 

X 

10-' 

240 

8.890 

X 

10-' 

250 

1.291 

X 

1 r\- 

260 

1.822 

X 

10-' 

270 

2.506 

X 

10-’ 

280 

3.370 

X 

10-' 

290 

4.440 

X 

10-' 

300 

5.742 

X 

10-- 
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APPENDIX  C 
THEORETICAL  REVIEW 


V i b ra_t_i  Ro ta  to^r  “'ode  1 

The  total  energy  of  the  molecular  vibrating  rotator  mooel  [«>-"-• 
1950]  is  the  sum  of  three  component  parts;  the  electronic  energy  , 
the  vibrational  energy  and  the  rotational  energy  That  is, 

= A- , + A],  + a],  (C.l) 

or,  if  we  write  in  the  equation  in  wavenumber  units  (or  term  values), 

T = /’  +;  + /■■  (C.2) 


where 


r r 

P = Planck's  constant,  6.625  x 10~‘^^  erg-sec. 

• - The  velocity  of  light,  2.99793  x 10‘°  cm/sec. 

The  wavenumber  of  the  s()0( tral  lines  corresponding  to  the  rotation 
between  iwu  elecfr'onic  states  ar'e  <jivf'n  I)y 

■ - M ”)  ^ (.  ■')  (■-  ' 

where  th<’  s i ngl  e- pr  ined  li'fters  relet'  to  the  ut>|>er  state  and  t tie  di",;  It 
primefi  nq.-r-  ’'i  the  lower  st.itf'. 

ih  ftp.  iMp'c,  , , iileeat  ion  is  n.ule  of  t h(>  possitile  transit'"'  . 

tn'twen  * he  r. 'at  ion. il  .t.ites  for  wtiuti  the  vitaation.  1 states  ah'* 
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electronic  states  are  given.  That  is,  only  the  rotational  transitions  of  0. 
from  v-O)  to  v=l).  The  first  and  second  terms  of  Equation  (C.3) 

can  be  combined  into  a single  term,  . , which  is  a constant  for  a specific 
vibrational  transition  and  is  called  the  band  origin.  Equation  (C.3)  can 
be  rewritten  as; 


r 


where 


' = the  total  angular  momentum  of  upper  level. 

. = the  rotational  term  of  a state  for  a given  >1  ’ . 

b , . = molecular  constants,  which  are  listed  in  Table  C.l. 

l’  0 

In  the  state,  which  means  A = 0,  .9  = 1,  .A  may  be  only  very  weakly 

coupled  to  the  internuclear  axis.  This  weak  coupling  of  S to  the  inter- 
nuclear  axis  is  the  characteristic  of  Hund's  case  {b) . If  we  designate 
the  total  angular  momentum  apart  from  spin  to  k,  k - .V  + A.  The 
quantum  number  of  can  have  the  integral  values; 


A = A , A + 1 , A + 2 , . 


(C.7) 


The  quantum  number  of  can  have  the  integral  values; 


= (a  + {r  + .■  - 1),  (a  + S - 2),....,  \k-s\.  (C.8) 

In  the  present  case  for  state,  A = 0,  5 = 1,  and  Equations  (C.7) 
and  (C.8)  become 


A'  = 0,  1 , 2,  3, 

J = (1),  (0,1,2),  (1,2,3),, 


(C.9) 

(C.IO) 


where  each  group  of  ,/'s  corresponds  to  A'  = 0,1, 2, 3, The  vector 

diagram  of  Hund's  i_ase  {b)  is  given  in  Figure  C.l.  In  this  figure 
the  nutation  of  the  figure  axis,  represented  by  the  broken-line 
ellipse,  is  much  faster  than  the  precessions  of  a and  and  r , re- 
presentf’d  hy  the  solid-line  ellipse.  For  A = , ( v)  is  perpen- 

diLular  to  ti)e  internucli-ar  axis. 

Tht>  rotational  ten"  of  c.ich  state  for  a given  (luantum  number  ■ 
ari>  qi  V'-n  t)y  ■>  , [ 1 936] : 

•1)  • (.  *3;  . - - y(2A+3)'^  , ( ■ + 1 ) 


1 


(c.n) 


:.•,(  .•)  = ^A-(^-+l)  - (2A-1)  P_^  \ +v{2A--l  )•/',/+>■ -2XJ- 

where 

/■'i(a)  = the  rotational  term  with  = A+l. 

A'2(A')  = the  rotational  term  with  ,/  = A'. 

/•’3(A)  = the  rotational  term  with  = A-1. 

B . D = molecular  constants  which  are  listed  in  Fable  C.l. 

o’  y 

X,  y = splittinq  constants  which  are  listed  in  Table  C.l. 

Transitions  of 

From  quantum  mechanical  theory,  the  transition  ♦ A'*2-) 

represents  a maqnotic  dipole  transtition.  The  selection  rules  that  *5 

hold  for  this  type  of  transition  are  listed  as  follows: 

(■■)  For  the  quantum  number  , of  the  total  angular  momentum  . 

= 0,  • 1,  with  the  restriction  .'  = 0 -'  = 0. 

(■■■■)  For  *he  quantum  number  ■'  of  the  total  angular  momentum  apart 
from  spin, 

0.  • 1 

and 

1.  forbidden  for  - ' transitions. 

(■'■)  lor  maqni-tii  di()ole  radiation,  the  rules  for  the  symmetry 
proper'* y are: 

t..*  «..»  -.fl 

and 

From  abov('  selection  rules,  four  tiranches  are  possible  for  each 
oiven  (juantum  numtier  of  upper  '.tate,  . ' . They  are: 


it 


1 

1 

f 
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; Table  C. 1 . Rotational 

[.'w/A-o.'f;', 

constants  of  the  a '):"  and 
1948]  •' 

■ *):'*’  states  of  0 . . 
9 

Constant 

1 

state 

(! 

(cm-*) 

/ state 

(cm- * ) 

1 

^ Vo  0 

13120.9080 

i 

! 1 

1566.3856 

i Ih 

1 .43777o 

1.39132  s 

: • 

1 .42197s 

1.37306s 

* Do 

4.91 3 X 10 

5.40,  X 10 

Lh 

4.82s  X 10 

5.45  s X 10 

A 

1.984 

T' 

-0.00837 

The  last  decimal  is 
within  one  unit  of  the 
less,  the  last  decimal 

printed  in  full  size  when 
next-to-the-1 ast  decimal, 
is  a subscript. 

the  accuracy  is  wel 1 
If  the  accuracy  is 

4 


(2)  Pq(-  '.A') 

(3)  >^Q(-/’,A’) 

(4)  RR(/^^') 


/■'(-' ',A') 


- A,(.  ',  A"  + l) 

- ..'-D 


= ')  - r , (.  '-1  , A'-l  ) 


where  , A'  refer  to  the  quuntufii  numbers  of  the  upper  level. 

The  superscripts  represent  the  changes  of  quantum  number  a when 
the  capital  letters  represent  the  changes  of  quantum  number  J.  P re- 
presents the  quantum  number  increased  by  1 . Q represents  equal  and 
R represents  the  quantum  number  decreased  by  1 when  the  transition  takes 
place  from  one  state  to  another. 

By  substituting  Equation  (C.ll)  into  (C.12),  the  vibrational  term 

values  of  the  lower  state,  .v’.T'  are  obtained.  From  Equation 

3 

(C.4)  and  (C.6)  the  wavelength  of  each  transition  which  corresponds 
each  line  of  transition  from  to  shown  in  Figure  C.2  can  be 
obtained.  In  Figure  C.2,  the  states  represented  by  broken  circles  are 
absent  because  of  zero  nuclear  spin.  The  wavenumbers  corresponding 
to  the  lines  in  Figure  C.2  are  generated  by  computer  and  are  listed 
in  Table  C.2. 
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Table  C.2. 

Wavenumbers 

(cm  ^ ) of  1 i nes 

in  (0-1 ) band 

of  02(*2-='Z). 

K‘ 

Branch 

'^Q 

Branch 

’'q 

Branch 

branch 

0 

11561.678 

2 

11570.026 

11571.906 

11557.887 

11555.807 

4 

11575.283 

11577.234 

1 1551 .701 

11549.691 

6 

11580.293 

11582.281 

11545.304 

11543.333 

8 

11585.054 

11587.069 

11538.675 

11536.731 

10 

11589.564 

11591 .603 

11531.806 

11529.885 

12 

11593.821 

11595.881 

11524.697 

11522.797 

14 

11597.823 

11599.903 

11517.345 

11515.46, 

16 

11601.567 

11603.667 

11509.750 

11507.890 

18 

11605.051 

11607.169 

11501 .91 1 

11500.070 

20 

11608.273 

11610.410 

11493.827 

11492.004 

22 

11611.228 

11613.383 

11485.495 

11483.690 

24 

1 1613.914 

1 1616.087 

11476.916 

11475.129 

26 

11616. 3'’rt 

11618.518 

11468.086 

11466.317 

28 

11618.464 

11620.673 

11459.006 

11457.253 

30 

11620.321 

1 1622.547 

11449.670 

11447.936 
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8 

7 

7 

RR 

1 15U'_..0G380 

0.531  732E  ♦ 1 I 

l*i 

8 

8 

7 

0 

KQ 

lir>8<^.069  18 

0.62/0121  111 

1 o 

8 

a 

9 

8 

F-0 

1 1538.67468 

0.656/151  *^1  1 

1 7 

B 

8 

9 

9 

FP 

11536. 73082 

0.7501521 +1  1 

la 

to 

10 

9 

9 

RK 

I 1589.56  563 

0.  46  537  IF  -M  1 

19 

10 

10 

9 

10 

KU 

11591 .60286 

0.528007L  ♦ 1 1 

?o 

10 

10 

1 

10 

FO 

11531  .8059  V 

0.5452  58E  + 1 ! 

?i 

10 

10 

1 

1 

F'P 

1 1529.88523 

0 . 608.5.141  ♦ 1 1 

27 

1 

12 

1 

1 

KK 

1 1593.82089 

0. 35  557  IF  + 1 1 

2S 

12 

12 

1 

12 

Kll 

1 1595 . 88098 

0 . 3939591  f 1 1 

?4 

12 

1 2 

13 

12 

f 8 

1 1524 . 69664 

0 . 4025  541  ♦ 1 1 

12 

1 2 

13 

1 5 

F'l 

1 1522. 79/49 

0.441  '80t  1 1 1 

26 

14 

14 

1 3 

1 3 

K h 

1 1597.82260 

0. 240.’86l  * \ 1 

1 4 

1 4 

1 3 

1 4 

f U 

i 15VV . 90280 

0.  .'635  521  ♦ 1 1 

2 a 

14 

1 4 

l*j 

1 4 

P8 

1 151  7.34516 

0 . 2669  961  f 1 1 

09 

1 4 

14 

r. 

r. 

f k 

11515.46612 

0.  .'HV4H1E  F 1 1 

iO 

16 

16 

ir> 

r. 

KK 

11601.  ..6652 

0.  1 4685  31  V 1 1 

i 1 

1 6 

16 

1 2 

K'i 

I 1605.666  8 4 

0 . 159  t 7.T  f I 1 

50 

1 6 

16 

1 7 

16 

f'l 

1 1509.  7500  5 

0 . 1 6<M  23F  1 1 

3 

\ 

16 

1 ? 

1 ■ 

K5 

1 1507.890  5 

0. 1 719891  + 1 1 

34 

IB 

IS 

1 7 

1 

KK 

1 I0O5.OM  r. 

0.81  IVV<'L+  10 

TS 

18 

le 

1 2 

18 

h'l 

) 160/ . 1 69  '5 

0.  .'23  721  MO 

3<S 

n: 

18 

19 

18 

^ U 

11501.91! 23 

0.8  M 101  ♦lO 

37 

li. 

10 

1 V 

19 

F F 

1 1 500 .0  /O.  4 

0. 92959  If  MO 

SU 

20 

20 

1 9 

19 

F.  K 

1 1608. .;2  V 

0 . 40  / H2  11  MO 

59 

20 

20 

I 9 

20 

116  10. 4096  1 

0 . 4 54fiv  11  ♦ 1 0 

40 

20 

20 

' 

2t< 

f 11 

1 1 493 . U2665 

0. 4 5.:.^4al  F 10 

41 

20 

20 

2 1 

2 1 

1 F 

1 1 49?  00  550 

0 . 45.1  1 6M  ♦ 10 

4? 

22 

21 

2 1 

KK 

1 I6l  1 . 22  769 

0. 1 H657M  F 1 ' 

4 5 

22 

2 2 

. 1 

•>  -» 

F>8 

1 161  3. 38  5 1 4 

0.  197 79  If  ♦ U' 

4 4 

2 2 

22 

2 3 

» 1 

f tJ 

1 1 495  4 t 

0. 195o28l  M • 

4C 

22 

22 

23 

5 

F F 

1 1 483  69(>  1 4 

0..’(5'..  .’9k  ♦ 10 

4a 

24 

2 4 

. 5 

2 3 

F>K 

1 1 6 1 3 . M 4 f 

0 . 779. 81.  k ♦ 09 

4 ■’ 

2 4 

. 4 

. 5 

. 4 

\i. 

1 1 616.086V  , 

W.021H85F  K)/ 

4 3 

24 

2 4 

2 4 

F a 

1 1 4 /6  . V f f,0  / 

0 . 80V  4 11  « ' 

4 •• 

.'4 

. 4 

2’, 

( F 

11  4 75.  1 . 

0 . .54  >■  .041  F . ' 

V. 

26 

2', 

Kk 

116  16.3.  \6  1 

’ 1 . / 1FI6I  ♦ ' 

•, 

.6 

. 

26 

KiJ 

1 16  13.51808 

. • 1 . 4 5 1 t 

. A 

> ; 

2A 

F U 

1 i 4^8.00.',  < 

”.1 

26 

.'6 

2/ 

1 » 

1 1 466  . 5 1 • , 

’ . 5 . . 1 5 /,  4 • F 0 ■ 

; 

20 

2s 

2 / 

2 / 

F'  > 

1 1 6 1 H . 4^  4 4 .» 

. 1 ..  tlUiVF  . . • 

'28 

2- 

. 8 

r.  ) 

1 I 62'»  . . 

'M  F 1 . 

28 

28 

29 

2H 

F 1 

1 1 -v  v.  • 

f . 1 >6  15(1  ♦ ■ ^ 

r'H 

29 

2 V 

r t 

114',  . . 5 . .1 

. ' 1 1 . . 4F  4 . ■ 

■ 4 

to 

29 

‘W 

F 

1 1 /.2'  .5.1 

M 5 1 ■ • ' 

, i 

( . 

29 

mF 

1 . ‘..5’. 54  .*•  . 

A i.  .i  *. 

5i» 

t) 

n 

5«' 

r M 

11449.6  .4 

. > . M 1 «1 1 4 , . 

!Ht  fOMUUlNii  PA1»'^  IS  5TUM  I.  IN 
THf  Tf  M>  I K.MURk  - IVti.OO 


NN 

0 Of  IRAN'S. 

fvTl  . 1 NT  1 NSl  T Y 

i 

1 r.M  .678  ^A 

70.47694 

■J 

1 l*.70.0?59fj 

V. 64774 

3 

1 15/1 ,905  76 

71 . 71800 

A 

1 155  88674 

76.44/87 

e 

1 15!, 5. 80654 

3U  . 440B6 

6 

1 1!. 75. 7831  1 

75. 10783 

/ 

1 I577.733f,7 

35.58037 

B 

11551. 70C53 

39.50371 

9 

1 1549.691 1 1 

49.8/339 

10 

1 1580.79788 

33.47356 

: 1 

1 1587.78069 

4 1 . 80096 

12 

1 1545 . 30470 

44 . 71404 

13 

11543. 33778 

57.96777 

1 4 

1 1585.05380 

34 . 4 3740 

ir» 

1 1587.069)8 

40.60777 

16 

1 I538.674.s8 

47. 57564 

17 

11536.73087 

48.5  '617 

IB 

1 1589,56363 

30.0O56H 

19 

1 1591 .60786 

34.  19118 

20 

1 1531  .80599 

35 . 30695 

71 

1 1579.8857  3 

39. 397/7 

">  -» 

1 1593.87089 

7. . 89557 

73 

11595.88090 

75 . 5108H 

74 

1 1574 . 69664 

76 . 0661 4 

7% 

1 r-77 . 79749 

78 . 60750 

76 

1 r.9  ■ . 07760 

1!,.  559  75 

7/ 

o 

> 

> 

1 7 • 06506 

7B 

1 1 !.  1 * . M 5 1 6 

1 . 78937 

J9 

11515. 466 1 7 

1 8 . 745  39 

30 

11601 .566' 7 

9 .50948 

31 

1 l6'<  3 . ^665  I 

10.30753 

37 

1 1 509 . 7500  I 

10.368  ^9 

33 

1 1507. 89».  <- 

1 1 . 1371R 

34 

1 1605.051 15 

5.75805 

3^ 

1 1607. 16935 

5 . 64  7 76 

3© 

1 1501  .9  M 78 

5 .64607 

3’ 

1 1500.07074 

6.0195.0 

3H 

1 1608 . 7 77  79 

.64085 

39 

1 1610 . 4.>96  1 

7.01M!, 

40 

1 1 49  I .0766*. 

7 . 79967 

41 

1149.'  .(>035.1 

7 .96686 

47 

11611. 77  "av 

1 . ?0H 1 / 

4 1 

1 16 1 3 . 38  U 4 

1 . . 8080 

44 

1 1 485 . 4954  1 

1 . 766 /V 

45 

1 1483 . 69014 

1 . 33576 

46 

11613.91437 

0 . f,O450 

47 

1 1616,08695 

0.53771 

48 

11476.91607 

0.57390 

49 

1 14  75. 17868 

0.55010 

50 

11616.37761 

0. 1975  ' 

51 

1 1618.51808 

0 . 70737 

57 

1 1468.00639 

0 . 1 98  76 

53 

1 1 4w66. 31607 

0.7074 J 

54 

11618. 46444 

0.06777 

55 

1 1670.677/9 

0.0/043 

56 

1 1 459 . 00567 

0.068/7 

5/ 

1 t 457 . 75374 

0.0/160 

58 

M 670 . 371  1 7 

0.07153 

59 

1 1 677 . 5466 1 

0,07747 

60 

1 1 449 .67004 

0.07184 

61 

1144^ . V3S5  , 

0.077/1 

6? 

11671. 893 1 9 

0.00637 

6 3 

1 1674. 1358  7 

0.00657 

64 

1 1 440.07HVI 

0 . 006  36 

6*. 

1 1 4 38.  361  5 , 

0 . 00660 

66 

1 1 67  < . 176  1 H 

O.OOl 70 

67 

1 1675. 4 r,V4 

O.OOl 76 

6H 

1 1 430. 7.  9,  4 

0.(’0l  70 

69 

1 1 4. 'H. 5790. 

O.OOl 76 

70 

1 16. ’4  . 1 656.' 

0 . Ov'047 

71 

1 16. V,  . 44.'*. 7 

0 . (H^U4  3 

7 7 

114  0.11731 

0 . (.00  4 7 

1141  '.43473 

0.00043 

74 

1 1 ' '4  . H56  9 

O . (hK'  1 0 

7», 

M 1M-’  ' 

O.-'OOlO 

76 

1 1 40v  . ’4  ■>.  - 

O . OOOOV 

r ’ 

. 1 408 . 0 ’!.A9 

0 . ooo 1 0 

THf 

TOT  Ai  1 N 1 1 N‘,  1 T y - 

1000.00008 

D-5 


Program  SORT 


k FIF 

*rry;  soRr/r 

r FROGRAh  SORF  DATA  ]N  AN  INCRfc AGING  OKU* R 

mm NG ION  no;  100)  .ntrr ( 100'  f RiNT<  100) 
i:  ««4iInmit  IiaTm  ♦*♦** 

1 Al  l IllFFN^  'SiYG' , m-’^Gl  ' ) 

I.  I 

• Kl  OLK  ) . i04  )N  . II.R 
iJn.R)  N 

mi  I ( R ) n 

R I N I ( Is  ) -R 
1 K Ml 

n (Iv-  ’/).5G.  TGrUl 
I:  I-  t,  I 

M R 

,■■■)  I M M 

M I f I X ( F M / r . 0 / 

IF  ( M ) IHO.  T’’  T . I no 
I no  IK  R M 

I I 

MG  1 I 

I I’O  I I F M 

IF  (iiirrci  ) mi  I 1 .1) ' ivo.  ivo.i'40 

?4n  II  IiIF 'Ml  ) 

I R I N T ( I , 

I I Nil  I I ) 

Ml  F (I  Ml  F > J) 

I : I N I a ) - R I N I ( I ) 

NiUl  ) ■ NiK  J) 

MFI(J)  ri 

RtNI  (.))  12 
NU  < , I ) - IT 
I 1 h 

IF  (I  ) MO.  1 VO.  I 20 
I V.)  i I t I 

ir<  I 1 R ) 1 1 0.  M ().  ^0 
S I 1 CON  I I NIIF 

UR  M 1 ( I . ’.T.l ) 

FIIRNAICIX.  I F NFTKA  Illl.f  1V5  R’/IX.'THF  GAIA  ARE  MORIG  IN 
I *MV’,2* 

1 IX..  N .2X.  INIiFX’.4X.  UAOF  NO.  .GX.'RIN!  ' 

C ♦♦♦(MTF'in  GAIA  **» 

I, i'll  I OUF  FNi  Til  . G2.TG2') 

111  I ■ ' I \ , ■ 7 

UK  1 I F I 4 , M 4 . G I F F II  ' . R- 1 N I ( I ) 

7.T4  F DRNm  I ( .'AA  I 

UR  I II  < 1 . .M'.  1 .Nil , I > . GIF  F ( I ) .R  INT  ( I ) 

yy:,  i ui.m.t i i m . i 4 . f 1 4 . •• . f 1 4 . g > 

7 • ■ RON  I I NIIF 

I Ai  I on  oGi 

104  F Ol.flAI  ' I 4 . . F 14. 

r.Ti  I F < 1 I 
I NG 


A.'X.F  14.  A) 


Output  of  Program  SORT 


N 

INI»1  X 

UAUt  NO. 

RINT 

1 

7 ; 

1 1408.07569 

0.00107 

o 

1 1 409 . 740R9 

0.00104 

7 3 

1 1418. 43348 

0.00368 

4 

73 

1 1 400. 1 1657 

0 . 00356 

/.9 

1 1420.50900 

0.01171 

ft 

6 8 

11430.00904 

0.01130 

' 

6". 

1 1438.36155 

0.03471 

8 

64 

1 1440.07891 

0.03345 

V 

6i 

1 1447.93480 

0.09579 

10 

60 

1 1449.67004 

0.09009 

1 1 

«■,  V 

1 1457.05304 

0.04597 

1 ? 

S6 

11459.00560 

0.03583 

1.1 

T/3 

1 1466.31613 

0.58733 

1 4 

t ') 

11468.08564 

0.56137 

1 

49 

11475. 10868 

1 . 30330 

16 

48 

11476.91532 

1.04120 

1 7 

45 

11483.69014 

0.68535 

IB 

44 

11485.49467 

0.54669 

19 

41 

1 1490.00350 

5. 13204 

:>o 

40 

11493.80665 

4 . 84093 

?1 

3 7 

11500.07004 

9.08604 

'>') 

36 

11501.91108 

8.52044 

23 

33 

11507.89037 

14 .87556 

1’4 

32 

11509.75003 

13.84904 

• »<- 

29 

11515.46537 

00.46691 

r<i) 

28 

11517.34516 

00. 70180 

1’ 

25 

11500.79675 

31 . 19969 

:*8 

24 

1 1524.69664 

08 . 40806 

29 

21 

1 1529 .88503 

39.64413 

?0 

20 

11531 .80504 

35.53008 

31 

17 

1 1536.73080 

45.74563 

3? 

16 

1 1538.67393 

40.04766 

33 

13 

1 1543 . 33004 

47 . 33344 

34 

12 

1 1545 . 30400 

39.95806 

9 

11549.69036 

40.88787 

36 

8 

11551 .69978 

33.97060 

3 7 

5 

11555.80654 

30.06451 

38 

4 

11557.88674 

00. 19383 

39 

1 

11561 .67834 

17.00409 

40 

2 

11570.00598 

8.09594 

A\ 

3 

11571 .90576 

18. 00474 

42 

6 

11575.2B036 

0 1 .58680 

43 

7 

11577 . 03090 

30.59676 

44 

10 

11580.09013 

09.86848 

4f, 

1 1 

11580.07995 

37. 35483 

46 

1 4 

11585.05305 

30. 40601 

4 /■ 

15 

1 1587.06918 

18.01630 

48 

18 

I 1589.5636  1 

30 .19  ’18 

47 

19 

11591  .600116 

14 . 409  IV 

!-0 

nn 

1 1593.80015 

04.97013 

23 

1 1595.80098 

07.80048 

26 

1 1597.80060 

1 8 . 64883 

*^3 

27 

1 1599. 90080 

00. 45098 

^4 

30 

11601 .56650 

1 0 . 701 49 

Ijfi 

S I 

1 1 603 . 66684 

1 3 . 76  '4  0 

S6 

34 

1 1 605 . 0504 1 

'.93675 

31. 

11607. 16861 

8.50501 

''.n 

38 

1 1608.070/9 

•1 . 56808 

1/9 

37 

11610.40961 

4.87  10.0 

60 

47 

11611 .22769 

0 . 428§4 

61 

43 

11613.38239 

0.57485 

6? 

46 

1 1613.91360 

1 .1952* 

61 

47 

11616.08621 

1 . 26092 

f>A 

50 

11616.32761 

0.54526 

65 

54 

11618.46444 

0.23085 

66 

51 

1 1618.51808 

0.57284 

6 ’ 

58 

11620.32037 

0.09079 

68 

55 

11620.67204 

0.24167 

6V 

6? 

1 1621 .89319 

0.03319 

70 

59 

11622.54661 

0.094  .’6 

71 

66 

11623.17618 

0.01129 

72 

6l 

11624 . 13507 

0.03453 

71 

70 

11624. 16562 

0.00337 

74 

74 

11624.85629 

0.00105 

75 

67 

11625.43520 

0.01172 

76 

71 

11626.44252 

0.00370 

77 

75 

11627.15107 

0.84189 

i 

1 

] 

i 

1 


n n n n nr. 


Program  CNVLS 


♦rTY:  CNOLS/T 

C ///THIS  FROGRAM  CONVOLWt  THt  INfUTE  HATG  HY  A RESOLUTION-RES  CM-1 

.'//ODTrUT  IiATA  FORMAT  ( I . SUM  ) - - ( I . I NT  b N' . I I Y ) ■ - ( All . AA  ) 

///INfUT  IiATA  IiIMF.N;;iON 
tUMENSION  UNN( /7) .RINTT ( /?) 

CALL  lOF'ENT  SYS' f I'i9f,0') 

K=  7 7 

DO  11  I^l.K 

REAIK  A f20'J)UNN(  I ) .K  INTT  ( I ) 

202  F0RMAI(2A6) 

11  CONTINUE 

///THE  UAUENUMBER  CORRESKIND  10  ORIGIN 
U- 11400.0 

///OUTTUT  DATA  FILl 
CALL  OOF'ENl  'SYS'  . ' DIVST  ' > 

////RT SOLUTION-RES  CM  1 
////INTERUAL  OF  ONE  'IIDF  R1C..4.0 
RES  3.0 
HINT=RES*4.0 
lB-1 

C ///THE  UF'PER  NUMBER  OF  DO  LOOP  IS  THE  NUMBER  UK  OUTPUT  POINTS 

DO  110  I=1»SOO 
U-UTO.S 
SUM-0.0 
J-JB 
N = 0 

3 DUN-UNN<J)-U 

IF<DWN>S<-30f15 
5 IF(DUN+DINT)10f2S.25 
10  J^J+1 

GO  TO  3 

25  N=N+1 

IF (N- 1 )30,26f30 

26  JB-J 

GO  TO  30 

15  IF(DUN-D1NT)30.30>60 

30  X=nWN*3. 1416/RES 
RINT^RINTTC J) 

AX=ABS(X) 

IF<AX-0. 000001 )40f40.50 
?0  RAT  I 0=1.0 
GO  TO  55 

50  RATIO  STN(X)/X 

55  SUM  = SUMFRINT»RATIO*( 1 .0-( DUN/DINT  1 ♦»2)**2 
TF(  J-M  IO.6O16O 
60  IF (SUM ) 70f 80.80 
70  SUM=0.0 

80  URI TE( 4.320) I .SUM 
320  FQRMAT(A2.A6) 

no  CONTINUE 

CALL  OCLOSE 
CALL  EXIT 
END 


Program  DTPRS 


tSYSIDTPRS  MT6S\;\:tiTPRS/A 
*TTY: < DTPRS/T 

[ilhENSIUN  UINN(  70)  . R I NT  T ( 70  > , R ( 70  > 

RFAFK  1 f 70?)N,(  l.NAMt 

202  f GRMAT(  RLOR  IlLE  NO.  = ' . 1 2/ ' OUTPUT  FILE  NAME  =• 
RC=  2047+ 1 

CALL  tUSRLOf  i'.0.0.20f0.  RC.IF.R) 

URITEL 1 p910) 

910  EORMATL 'CHfl  R ' ) 


i I AH. lUAL) 


EORMATL 'CHfl  R ' ) 

IAD=^-280 
DO  11  1 - 1 . 70 

CALL  DATAHLO.:’.  IADf  IVAl  ) 

XP  = K>AL 
lAD  lAI'tl 

CALL  liATASLO.2.  lAD.  lUAL) 

XMAN  lOAL 

RINTI  ( I ) =-XMAN»  L2.0»*XP) 

iad-tad  + ;l 

CONTINUE 

URITf  L 1 f340) ( I pRINTTL I > fI  = 1 ,6) 

F0RMAr(2(4Xf I2.E16.6) ) 

///INSTRUMENTAL  RESPONSE  CORRECTION  AND  NORMAL  IZATION 
///SET  THE  START  DATA  POINT 
ISTR-2V 
DMINI =0.0 

///'DAIR'  IS  THE  FILE  FOR  INSTRUMENTAL  RESPONSE 
CALL  lOPENL  SYS' . 'DAIR' ) 

READL  4. 302) (R( I ) f I 1 f 40 ) 

F0RMAT(Eri.3) 

DO  44  I 4L  f70 

R ( I ) = i . 00 

CONTINUE 

TOTL  0.0 

DO  33  I ISTRf  70 

RINTI < I )=RINTT< I )-DMINI 

RINTTL I )^RTNTT( I ) 7R ( I ) 

TOTL  - TOTL+RINTTL I ) 

20NT 1 NUE 
IEND-70-ISTR+1 
DO  22  I=1fIEND 
N I+ISTR-1 

RINTTLD-LRINT  T ( N )»  1 000 . 0 ) /TOTL 
FN  = N 

UNNL I )=1I457.46+FN*2.23 
CONTINUE 

WRI TEL  1 f360) L I .UNNL  DfRI NT TLI). 1 = 1.6) 

EORMATL I3f3XfF10.2f3XfE12.4) 

///f'ONOOLVE  THE  DATA  PY  A RESOLUTION=RES  CM-1 
URITEL 1 .910) 

///SET  R=THE  NO.  OF  LINES  OF  LINE  SPECTRUM 
R = 42 

///THE  UAOENUMPER  CORRESPOND  TO  ORIGIN 
U=1 1400.0 

///OUTPUT  DATA  FILE 
CALL  OOPENL 'SYS' fFLNAME) 

/.  VRESOLUTION=RES  CM-1 
/.//INTERNAL  OF  ONE  SIDE=  RES*4.0 
RES  3.0 
DINT=RES*4 . 0 
JP  - 1 

///THE  UPf ER  NUMBER  OF  DO  LOOP  IS  THE  NUMBER  OF  OUTPUT  POINTS 
DO  110  I 1.500 

U-U+C.5 
SUM  =0.0 
J«JB 

NO 


10 


3 DUN-UNN(J)-W 
IF(DUN)5t30rl5 
5 IF  (riUN+niNT)10,23»25 

10  J J+1 

(ill  ro  3 

25  N^N+1 

IF(N  1)30.26.30 

26  JF  J 

r,n  10  30 

15  IF  (liUN-rHNT)30,30.60 
30  X-DUN)*3.  M16/Ki;S 

FINT-KINITI J> 

OX  APS(X) 

IF  (AX-0.000001  )40.'10.50 
^0  RATIO-- 1.0 
on  ru  55 

50  RATIO- S1N(X)/X 

55  3un^  SOMFRINI*RATIO»(  1 . 0 ' ( IiUN/IiINT  ) **2  ) ** 
IF ( J R ) 10.60.40 
60  IF(SUM)70.80.90 

70  SUM ^0.0 
BO  UK  I TF  ( 4 .320) I .SUM 
320  FORMAT (A2.A6) 

110  CONMMl-K 

I All  OCl  OSE 
I Al  L EXIT 
1 Nil 
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APPENDIX  E 

VALUES  OF  INSTRUMENTAL  RESPONSE 

The  values  in  the  following  table  are  stored  in  data  filed  'DAIR  DA.' 
Total  number  of  values  is  70.  The  rest  of  them  are  equal  to  1. 


No. 

Relati ve 
Response 

No. 

Relative 

Response 

1 

0.715 

21 

0.904 

2 

0.725 

22 

0.912 

3 

0.735 

23 

0.920 

4 

0.748 

24 

0.928 

5 

0.769 

25 

0.935 

6 

0.771 

26 

0.942 

7 

0.781 

27 

0.949 

8 

0.790 

28 

0.9t>u 

9 

0.801 

29 

0.962 

10 

0.810 

30 

0.969 

11 

0.820 

31 

0.974 

12 

0.830 

32 

0.978 

13 

0.840 

33 

0.982 

14 

0.849 

34 

0.985 

15 

0.857 

35 

0.988 

16 

0.868 

36 

0.991 

17 

0.874 

37 

0.994 

18 

0.881 

38 

0.996 

19 

0.889 

39 

0.998 

20 

0.896 

40 

0.999 

APPENDIX  F 
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Figure  F.l.  The  observation  of  zenith 
eiaission  of  02(.'‘?:t)  of  atmospheric  (0,1) 
band  from  rocket;  the  points  refer  to 
spectrometer,  and  the  dash  line  refers  to 
photometer.  The  solid  line  is  the  adjust- 
ed curve  to  make  the  intensities  agree. 
[Wallui-'e  and  Hunten,  1968]. 

OVERHEAD  EMISSiOH  (HR)  for  SI 
I 10  100 


OVERHEAD  EMISSION  (liR)  for  PI 


Figure  f.2.  Zenith  emission  rates  observ(>d 
by  photometer  P.  and  spectrometer  SI. 

• and  ii..  ' . >:,  1 968] . 


(11,564) 


DISTRIBUTION  LIST 


I Department  of  Defense 

' Director 

Defense  Advanced  Rsch  Proj  Agency 
Architect  Building 
^ 1400  Wilson  Blvd. 

■ Arlington,  VA  22209 

I 

OICY  Attn:  LTC  W.  A.  Whitaker 

j OICY  Attn:  STO  Capt  J.  Justice 

i OICY  Attn:  Major  Gregory  Canavan 

! Defense  Documentation  Center 

I Cameron  Station 

i Alexandria,  Va  22314 

(12  copies  if  open  publication, 
otherwi se  02  CY ) 

12CY  Attn:  TC 


Director 

Defense  Nuclear  Agency 
Washington,  D.C.  20305 


03CY 

Attn: 

STTL 

Tech  Library 

OICY 

Attn: 

STS  I 

Archives 

03CY 

Attn: 

RAAE 

Charles  A.  Blank 

OICY 

Attn: 

RAAE 

Harold  C.  Fitz,  Jr 

OICY 

Attn: 

RAAE 

Maj  John  Clark 

Dir,  of 

Defense 

Rsch 

& Engineering 

Department  of  Defense 
Washington  D.C.  20301 

OICY  Attn:  DD/S^SS  Daniel  Brockway 

OICY  Attn:  DD/SfSS  Richard  S. 

Ruffine 

Commander 
Field  Command 
Defense  Nuclear  Agency 
Kirtland  AFB,  NM  8411b 

OICY  Attn:  FCPR 

Chief 

Livermore  Division  FLD  Command  DNA 
Lawrence  Livermore  Laboratory 
P.  0.  Box  808 
Livermore,  CA  94550 


Department  of  Army 

Commander/Director 
Atmospheric  Sciences  Laboratory 
U.S.  Army  Electronics  Command 
White  Sands  Missile  Range,  NM  88002 

OICY  Attn:  DRSEL-BL-SY-S 

F,  E.  Niles 

OICY  Attn:  E.  Butterfield 

DRSEL-BL-SY-R 

Commander 

Harry  Diamond  Laboratories 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

(CNWDI-Inner  Envelope:  Attn: 

DRXDO-RBH) 

02CY  Attn:  DRXDO-NP 

Commander  ^ .j 

U.S.  Army  Nuclear  Agency 
Fort  Bliss,  TX  79916 

OICY  Attn:  ATCA-NAW  J.  Berberet 

D^aj_tjTient_of  f^vy 

Chief  of  Naval  Research 
Navy  Department 
Arlington,  VA  22217 

OICY  Attn:  Code  427  CDR 

Ronald  J.  Oberle 

Commander 

Naval  Electronics  Laboratory  Center 
San  Diego,  CA  92152 

OICY  Attn:  Code  2200  1 Verne  L. 

Hi  1 derbrand 

OICY  Attn:  Code  2200  Ilan 

Rothmul 1 er 


OICY  Attn:  FCPRL 


H-2 


Director 

Naval  Research  Laboratory 


Washington,  D. 

C.  20375 

OICY 

Attn: 

Douglas  P.  McNutt 

OICY 

Attn: 

Code  7127 

Charles  Y.  Johnson 

OICY 

Attn: 

Code  2027  Tech  Lib 

OICY 

Attn: 

Code  7700 
Timothy  P.  Coffey 

OICY 

Attn: 

Code  7701 
Jack  D.  Brown 

OICY 

Attn: 

7750 

Darrell  F.  Strobel 

OICY 

Attn: 

Code  7750 
Paul  Juluenne 

Commander 

Naval  Surface  Weapons  Center 
White  Oak,  Silver  Spring,  MO  20910 

01 CY  Attn:  Code  WA501 

Navy  Nuc  Prgms  Off 


Department  of  the  Air  Force 

AF  Geophysics  Laboratory,  AFSC 
Hanscom  AFB,  MA  01731 

01 CY  Attn 
OICY  Attn 
01 CY  Attn 

AF  Weapons  Laboratory,  AFSC 
Kirtland  AFB,  NM  87117 

OICY  Attn:  DYT  Capt  David  W.  Goetz 

OICY  Attn:  SUL 

OICY  Attn:  DYT  LTC  Don  Mitchell 

Commander 

ASD 

WPAFB,  OH  45433 

OICY  Attn:  ASD-YH-EX 

LTC  Robert  Leverette 

SAMSO/SZ 

Post  Office  Box  92960 
Worldv;ay  Postal  Center 
Los  Angeles,  CA  90009 
(Space  Defense  Systems) 

OICY  Attn:  SZJ  Maj  Lawrence  Doan 


U.S.  Energy  Rsch  and  Dev  Admin 

Division  of  Military  Application 
U.S.  Energy  Rsch  & Dev  Admin 
Washington,  D.C.  20545 

OICY  Attn:  Doc  Con  for 

Maj.  D.A.  Haycock 

Los  Alamos  Scientific  Laboratory 

P.O.  Box  1663 

Los  Alamos,  NM  87545 

OICY  Attn:  Doc  Con  for 
R.A.  Jeffries 


Other  Government 

Department  of  Commerce 
Office  of  Telecommunications 
Institute  for  Telecom  Science 
Boulder,  CO  80302 

OICY  Attn:  William  F.  Utlaut 

OICY  Attn:  Glenn  Falcon 


Department  of  Defense  Contractor 

Aerodyne  Research,  Inc. 

Bedford  Research  Park 
Crosby  Drive 
Bedford,  MA  01730 

OICY  Attn:  F.  Bien 

OICY  Attn:  M.  Camac 


Aerospace  Corporation 

P.O.  Box  92957 

Los  Angeles,  CA  90009 


OICY 

Attn ; 

T.  Taylor 

OICY 

Attn : 

R.  Grove 

OICY 

Attn : 

R . D.  Rawcl i ffe 

OICY 

Attn ; 

Harris  Mayor 

Denver,  University  of 
Colorado  Seminary 
Denver  Research  Institute 
P.O.  Box  10127 
Denver,  CO  80210 

(Only  1 copy  of  Class  P.pts) 
OICY  Attn;  Sec  Officer  for 
Mr.  Van  ZYL 

OICY  Attn;  Sec  Officer  foi' 
Oavid  Mutaa'ay 


LKB  Kenneth  S.W.  Champion 
OP  John  S.  Garing 
OPR  Alva  T.  Stair 


H-3 


General  Electric  Company 
Tempo-Center  for  Advanced  Studies 
816  State  Street  (P.O.  Drawer  QQ) 
Santa  Barbara,  CA  93102 

05CY  Attn;  Dasiac  Art  Feryok 
01 CY  Attn:  Warren  S.  Knapp 

General  Research  Corporation 

P.O.  Box  3587 

Santa  Barbara,  CA  93105 

OICY  Attn;  John  Ise,  Jr. 

Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 

(All  Class  Attn:  Security  Officer) 

OICY  Attn:  T.  H.  Davis 

(Unci . only) 

03CY  Attn:  Neal  Brown 

(Unci,  only) 

Honeywell  Incorporated 
Radiation  Center 
2 Forbes  Road 
Lexington,  MA  02173 


Mission  Research  Corporation 

735  State  Street 

Santa  Barbara,  CA  93101 

OICY  Attn:  P.  Fischer 

OICY  Attn:  D.  Archer 

Photometries,  Inc. 

442  Marrett  Road 
Lexington,  MA  02173 

OICY  Attn:  Irving  L.  Kofsky 

Physical  Dynamics  Inc, 

P.O.  Box  1069 
Berkeley,  CA  94701 

OICY  Attn:  Joseph  B.  Workman 

Physical  Sciencer,  Inc. 

30  Commerce  Way 
Woburn,  MA  01801 

OICY  Attn;  Kurt  Wray 

R & D Associates 

P.O.  Box  9695 

Marina  Del  Rey,  CA  90291 


OICY  Attn:  W.  Williamson 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 

OICY  Attn;  Ernest  Bauer 
OICY  Attn:  Hans  Wolfhard 

Lockheed  Missiles  and  Space  Company 
3251  Hanover  Street 


Palo 

A1 

to,  CA 

94304 

01 

CY 

Attn: 

John 

B.  Cladis 

Dept. 

52-12 

01 

CY 

Attn: 

J.  B. 

Reagan 

D/52- 

12 

01 

CY 

At  tn: 

Billy 

M.  M'-f.orr;ac 

Dept. 

52-54 

01 

CY 

Attn; 

Tom  Janies 

01 

CY 

Attn: 

Roher 

t n.  Sears 

Dep. 

52-14 

01 

CY 

Attn; 

John 

f ai:  ('r' 

01 

CY 

At  Ln ; 

Mart  i 

n Wal  , ,'t.  5 

01 

CY 

Attn: 

R i . ha 

rd  0.  ..iiiinson 

Dept , 

52-12 

OICY  Attn:  f^orrest  Gilmore 

OICY  Attn:  Robert  E.  Lelevier 

R & D Associates 
1815  N.  Ft.  Myer  Drive 
11th  Floor 

Arlington,  VA  22209 

OICY  Attn;  Herbert  J.  Mitchell 

The  Rand  Corporation 
1 700  Main  Street 
Santa  Monica,  CA  90406 

OICY  Attn:  James  Oakley 

Science  Applications,  Inc. 

P.O.  Box  2351 
La  Jolla,  CA  92038 


OICY  Attn: 


Daniel  A.  Hamlin 


H-4 


Space  Data  Corporation 
1331  South  26th  Street 
Phoenix,  AZ  85034 

01 CY  Attn:  Edward  F.  Allen 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
' Menlo  Park,  CA  94025 

OICY  Attn:  Walter  G.  Chestnut 

OICY  Attn:  M.  Baron 

OICY  Attn:  Ray  L.  Leadabrand 

Stanford  Research  Institute 
1611  North  Kent  Street 
I Arlington,  VA  22209 

i 

OICY  Attn:  Warren  W.  Berning 

Technology  International  Corporation 
75  Wiggins  Avenue 
Bedford,  MA  01730 

OICY  Attn;  W.  P.  Boquist 

Utah  State  University 
Logan,  UT  84322 

OICY  Attn;  Doran  Baker 
OICY  Attn:  Kay  Baker 

' OICY  Attn:  D.  Burt 

OICY  Attn;  C.  Wyatt 

Vi  si  dyne,  Inc. 

19  Third  Avenue 

North  West  Industrial  Park 

Burlington,  MA  01803 

OICY  Attn:  William  Reidy 

OICY  Attn;  J.  W.  Carpenter 
OICY  Attn:  T.  C.  Degges 


